Modelling the interaction between induced pluripotent stem cells derived cardiomyocytes patches and the recipient hearts

Cardiovascular diseases are the main cause of death worldwide. The single biggest killer is represented by ischemic heart disease. Myocardial infarction causes the formation of non-conductive and non-contractile, scar-like tissue in the heart, which can hamper the heart's physiological function and cause pathologies ranging from arrhythmias to heart failure. The heart can not recover the tissue lost due to myocardial infarction due to the myocardium's limited ability to regenerate. The only available treatment is heart transpalant, which is limited by the number of donors and can elicit an adverse response from the recipients immune system. Recently, regenerative medicine has been proposed as an alternative approach to help post-myocardial infarction hearts recover their functionality. Among the various techniques, the application of cardiac patches of engineered heart tissue in combination with electroactive materials constitutes a promising technology. However, many challenges need to be faced in the development of this treatment. One of the main concerns is represented by the immature phenotype of the stem cells-derived cardiomyocytes used to fabricate the engineered heart tissue. Their electrophysiological differences with respect to the host myocardium may contribute to an increased arrhythmia risk. A large number of animal experiments are needed to optimize the patches' characteristics and to better understand the implications of the electrical interaction between patches and host myocardium. In this Thesis we leveraged cardiac computational modelling to simulate \emph{in silico} electrical propagation in scarred heart tissue in the presence of a patch of engineered heart tissue and conductive polymer engrafted at the epicardium. This work is composed by two studies. In the first study we designed a tissue model with simplified geometry and used machine learning and global sensitivity analysis techniques to identify engineered heart tissue patch design variables that are important for restoring physiological electrophysiology in the host myocardium. Additionally, we showed how engineered heart tissue properties could be tuned to restore physiological activation while reducing arrhythmic risk. In the second study we moved to more realistic geometries and we devised a way to manipulate ventricle meshes obtained from magnetic resonance images to apply \emph{in silico} engineered heart tissue epicardial patches. We then investigated how patches with different conduction velocity and action potential duration influence the host ventricle electrophysiology. Specifically, we showed that appropriately located patches can reduce the predisposition to anatomical isthmus mediated re-entry and that patches with a physiological action potential duration and higher conduction velocity were most effective in reducing this risk. We also demonstrated that patches with conduction velocity and action potential duration typical of immature stem cells-derived cardiomyocytes were associated with the onset of sustained functional re-entry in an ischemic cardiomyopathy model with a large transmural scar. Finally, we demonstrated that patches electrically coupled to host myocardium reduce the likelihood of propagation of focal ectopic impulses. This Thesis demonstrates how computational modelling can be successfully applied to the field of regenerative medicine and constitutes the first step towards the creation of patient-specific models for developing and testing patches for cardiac regeneration.Open Acces

Modelling the interaction between stem cells derived cardiomyocytes patches and host myocardium to aid non-arrhythmic engineered heart tissue design

Application of epicardial patches constructed from human-induced pluripotent stem cell- derived cardiomyocytes (hiPSC-CMs) has been proposed as a long-term therapy to treat scarred hearts post myocardial infarction (MI). Understanding electrical interaction between engineered heart tissue patches (EHT) and host myocardium represents a key step toward a successful patch engraftment. EHT retain different electrical properties with respect to the host heart tissue due to the hiPSC-CMs immature phenotype, which may lead to increased arrhythmia risk. We developed a modelling framework to examine the influence of patch design on electrical activation at the engraftment site. We performed an in silico investigation of different patch design approaches to restore pre-MI activation properties and evaluated the associated arrhythmic risk. We developed an in silico cardiac electrophysiology model of a transmural cross section of host myocardium. The model featured an infarct region, an epicardial patch spanning the infarct region and a bath region. The patch is modelled as a layer of hiPSC-CM, combined with a layer of conductive polymer (CP). Tissue and patch geometrical dimensions and conductivities were incorporated through 10 modifiable model parameters. We validated our model against 4 independent experimental studies and showed that it can qualitatively reproduce their findings. We performed a global sensitivity analysis (GSA) to isolate the most important parameters, showing that the stimulus propagation is mainly governed by the scar depth, radius and conductivity when the scar is not transmural, and by the EHT patch conductivity when the scar is transmural. We assessed the relevance of small animal studies to humans by comparing simulations of rat, rabbit and human myocardium. We found that stimulus propagation paths and GSA sensitivity indices are consistent across species. We explored which EHT design variables have the potential to restore physiological propagation. Simulations predict that increasing EHT conductivity from 0.28 to 1–1.1 S/m recovered physiological activation in rat, rabbit and human. Finally, we assessed arrhythmia risk related to increasing EHT conductivity and tested increasing the EHT Na(+) channel density as an alternative strategy to match healthy activation. Our results revealed a greater arrhythmia risk linked to increased EHT conductivity compared to increased Na(+) channel density. We demonstrated that our modeling framework could capture the interaction between host and EHT patches observed in in vitro experiments. We showed that large (patch and tissue dimensions) and small (cardiac myocyte electrophysiology) scale differences between small animals and humans do not alter EHT patch effect on infarcted tissue. Our model revealed that only when the scar is transmural do EHT properties impact activation times and isolated the EHT conductivity as the main parameter influencing propagation. We predicted that restoring physiological activation by tuning EHT conductivity is possible but may promote arrhythmic behavior. Finally, our model suggests that acting on hiPSC-CMs low action potential upstroke velocity and lack of I(K1) may restore pre-MI activation while not promoting arrhythmia